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O ■ Abstract 

^ ■ The permeation, rejection, and transport of electrolytes in water-filled nanopores are critical to 

. ion current gating and desalinalion processes in synthetic porous membranes and the functions of 

> . 

biological ion channels. While the effects of confinement pore polarizability, and discrete channel 

^ ' charge sites have been much studied, the potentially dramatic impact of dipole-lined synthetic pores 

■ on electrolytes has seldom been addressed. Dipole layers naturally occur on the interior surfaces 

O . of certain nanopores, leading to intrinsic preference for cations or anions. This preference can be 

exploited when the membrane surface is functionalized differently from the pore interior or when 

^ ' there are alternating dipolar /nondipolar stretches inside a long pore. The dipole- ion interaction is 
l> 

lO . asymptotically unscreened by water, leading to ionic, charge segregated, insulating behavior that 
(N 

'nI" I can block ion transport, and potentially novel current-voltage (I — V) characteristics. 
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Water-filled nanoporcs exhibit a remarkable ability to exclude or selectively transmit 
ions.^"^ They have important biological functions, and are potentially useful for desalina- 
tion and other applications. The effects of dielectric mismatch/ nanopore diameter,^ net 
charge/ and polarizability^'^'^ on ion permeation and rejection have been much studied. 
Surprisingly, while the celebrated KcsA channel selectivity filter relies on carbonyl group 
dipoles pointed radially into its narrow passageway for its function,^ dipolar decorations in 
synthetic nanopores have received less attention. A recent force field-based model has how- 
ever shown that synthetic silica nanopores contain dipoles on its interior surface and exhibit 
an intrinsic preference for cations — even at the pH of zero charge.^ This work elucidates 
the dramatic effect of dipolar layers in nanopores on electrolyte permeation, which may lead 
to novel current-voltage (/ — V) characteristics, and discusses pore geometries where this 
effect may be realized. 

An infinitely long cylindrical pore of radius Rd with radially aligned surface dipoles of 
magnitude d and uniform surface density rid exerts a uniform electrostatic potential 0d = 
—ATrrjdd in the pore interior, independent of Rd or radial position of ion R. This intrinsic 
0d|e| is predicted to be ~ —24 kcal/mol in silica pores,^ and exceeds -100 kcal/mol inside 
the KcsA filter. Here e is the electronic charge. A uniform dipole layer thus induces no 
electric field and does not interact with uncharged H2O molecules. Hence its interaction with 
ions inside the pore is unmediated (unscreened) by the confined water. This is consistent 
with our hitherto puzzling observation that the free energy preference for Na+ over Cl~ in 
the infinitely long silica pore of Ref. 6 is fairly independent of the presence of water. Thus 
the simple uniform dipole density in our model apparently captures essential nanopore-ion 
interactions in nanoporous silica, where H2O molecules do form hydrogen bonds with surface 
silanol groups. This model should be pertinent to other inorganic nanopores too. 

We next use two models to demonstrate how this intrinsic (pd niay be exploited: (1) 
an infinitely long "tube" geometry, with repeating unit cells that have alternating L/2- 
long stretches where dipoles surface densities exist or vanish (Fig. la); (2) a "membrane" 
geometry with dipoles in the pore interior only (Fig. lb). Regions with low local dipole 
density can potentially be achieved by functionalizing silica pores with amine groups.^ Near 
the junction, z = Zo, between dipolar (0 = —A-nrjdd) and non dipolar (0 = 0) regions, the 
pore-axis electrostatic potential is 

-2^r)ddl ^"'""^^^'.^ + 1] (1) 



0(2;, R) is fairly independent of R inside the pore but vanishes rapidly for R > R^. An 
electric field parallel to the pore interior surface materializes/^ localized within a distance 
Rd of the entrance. This field is locally screened by water molecules. 

Note that in the membrane model, if the membrane surface has the same dipole density 
as the pore interior surface, it exerts an equal and opposite (pdi^, R — 0) that cancels Eq. l}^ 
Hence, to exploit dipolar effects, the two surfaces must be functionalized differently. 

To examine ion permeation into water-filled dipolar pores, we apply the atomistic H2O, 
Na"*", and Cl~ force fields used in Ref. 6 to compute the potential of mean force W{z). Pore- 
water/ion interactions are described by analytical, least square fit approximations to dipolar 
electrostatic and Lennard- Jones terms. ^° The pores are hydrophilic, with ~ 12 A diameter 
and 2.7 H2O density. Infinitely long but laterally isolated tubes (Fig. la) are simulated 
using at least 10^ passes of canonical ensemble Monte Carlo (MC) runs. Membrane models 
(Fig. lb), periodic in all 3 directions with lateral unit cell lengths of 40 A, are simulated for 
at least 200 ps using molecular dynamics (MD) /particle mesh Ewald methods. Simulation 
cells contain two unit cells doubled in the z direction. Because Na"*" orders water in nanpores 
in a way incompatible with the periodic boundary conditions used,^^ our simulation cells 
contain two unit cells doubled in the z direction for convergence. 

Figure 2b depicts W{z) for a 0d = —72 kcal/mol and L/2 = 30 A tube. The membrane 
geometry leads to qualitatively similar results (Fig. 2c). These geometries yield AWi = 
W{L/2) - W{0) = -8.1 and -12.5 kcal/mol. Compared with AWi = 0d=-72 kcal/mol in 
the absence of water, they exhibit effective dielectric constants of 8.9 and 6.0, respectively. 

However, unlike transmembrane ion channels such as KcsA,^ material pores can have 
lengths on the order of 1 //m. Table I shows that Al^i in the tube geometry steadily 
increases with L/2 and is roughly proportional to (f)^', at L/2 = 120 A, e ~ 2.3. This is 
readily explained as follows. Unless the localized electric field induces a ferroelectric phase 
transition (not the case here), AWi is predominantly screened by H2O near the junctions 
where the electric field at its maximum. For large L, these localized H2O resemble a point 
dipole whose interaction with ions vanishes as l/\z — Zo\'^- Thus, the attraction of cations 
deep inside dipolar regions is asymptotically unscreened by water. AWi analysis for the 
membrane geometry is complicated by its non-zero value even when the dipoles are absent in 
membrane slabs, and by the 0(1) kcal/mol statistical noise, but AWi follows a qualitatively 
similar trend there. 



screen the ion-surface dipole interaction. Here nanoscale pores with Rd < 10 A may exhibit 
unique properties. (1) We consider 1.0 M NaCl electrolytes in the membrane geometry. The 
hydration free energy of the Na"*" in the middle of the dipolar pore is reduced (last hues of 
Table I) because other Na+ ions now enter the nanopore, while Cl~ seldom do; the resulting 
electrostatic potential increase makes cations less favorable. ^'^ But unlike macroscopic pores, 
even the 1.0 M NaCl electrolytes do not yield complete screening in this membrane geometry. 
This is possibly due to confinement efi:ects. Results from 2.7 and 5.1 ns MD trajectories 
with different initial configurations are used here because ions diffuse rather slowly. (2) 
Azobenzene^^ or similar functional groups can be placed at the junctions where the axial 
electric field is largest, further restricting water density there and hindering screening. (3) 
Nanopores with alternating dipolar and non-dipolar surfaces may result in an ionic, charge 
segregated, insulating behavior. 

We use MC to demonstrate that such a "phase separation" occurs with a confined 0.5 M 
NaCl inside an infinitely long pore with unit cell L — 240 A and 0d|e| — —72 kcal/mol (see 
Fig. Ic). The effective dipole-ion screening for this L is 2.3 (Table I). This Na+ concentration 
is qualitatively consistent with our simulation of permeation of 1.0 M NaCl into the nanopore 
in the membrane geometry. ^° To accelerate ion sampling, water is here approximated as a 
dielectric with e„ = 10, likely slightly underestimated because this is obtained by considering 
a contact ion pair in the tube.^'^ We find that ion transport is blocked even upon applying a 
0.0025 V/A electric field screened by e„ = 10. We stress that this insulating behavior arises 
from (1) the almost- unscreened ion-pore dipole interaction dominating the better-screened 
{e^ — 10), ion-ion coulomb interactions, and (2) the alternating dipolar regions leading to 
a complete absence of cation/anion along lengths of the nanopore. Future simulations of 
I — V characteristics, taking into account leakage of Na+ through the nondipolar regions at 
long times, will be profitable. 

In conclusion, dipole layers naturally occur in certain nanopores, leading to intrinsic 
preference for cations or anions. This preference can be exploited when the membrane surface 
is functionalized differently from the pore interior, or when there are alternating dipolar/non- 
dipolar stretches inside a long pore. The dipole-ion interaction is asymptotically unscreened 
by water, potentially leading to novel ionic, charge segregated, insulating behavior that can 
block ion transport. Our finding may lead to new synthetic ion-exclusion membranes. 
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120 240 360 480 

FIG. 1: (a) Cylindyrical nanopore unit cell for the "tube" geometry. Ions and water are excluded 
from the crossed regions. Radial dipole moments ("+/-") exist on alternating stretches along the 
pore surface, (b) "Membrane" simulation unit cell. Blue and red/white spheres depict Na"*" and 
II2O respectively. Dipoles exist on interior pore surface only, (c) Charge segregated insulating 
behavior in tube geometry (see text). Na"^ and Cl~ are depicted as blue and green spheres. 
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FIG. 2: Na"*" potential of mean force W{z) kcal/mol along the pore axis, with (j)(>\e\ = 71.7 kcal/mol, 
L = 60 A, and a unit cell doubled in the ^-direction, (a) No water; solid, dotted, and dashed lines 
refer to ions at radii R/Rd=0, 0.7, and 1.8 respectively, (b) Water-filled tube geometry, (c) 
Water-filled membrane geometry. 
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TABLE I: Na"^ potential of mean force (APFi) at z = L/2, relative to its value at z = where 
PS 0. Energies and lengths are in units of kcal/mol and A. *Computed in 1.0 M NaCl electrolyte, 
not pure water. 



